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Abstract. The theoretical K factor, describing the difference between the leading and 
higher order cross sections, has no precise definition. The definition is sensitive to the 
order of the fit to the parton densities and the number of loops at which a.s is evaluated. 
We describe alternate ways to define the K factor and show how the definition affects 
its magnitude and shape for examples of hadroproduction of W+ bosons, Drell-Yan 
lepton pairs, and heavy quarks. We discuss which definition is appropriate under 
certain circumstances. 
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1. Introduction 

It has been clear for many years that only a leading order, LO, evaluation of quantum 
chromodynamics, QCD, cross sections is inadequate to describe Drell-Yan, heavy quark, 
and jet production. For example, at LO lepton pairs from the Drell-Yan process and heavy 
quark pairs are produced back-to-back with zero transverse momentum, pr, of the pair. 
Thus there is no nonzero Drell-Yan or QQ pair pr at LO. In addition, LO cross sections 
underestimate the measured total cross sections by up to a factor of two or more. 

Next-to-leading order, NLO, evaluations of these cross sections removed many of these 
inadequacies, especially concerning finite Drell-Yan or QQ pair pr since light quark or 
gluon emission at NLO keeps the pair from being perfectly back-to-back. In addition, 
processes such as gg — > gg* gQQ produce high pr QQ pairs since the pair takes the 
entire pr of the excited gluon. The magnitude of the Drell-Yan cross section at NLO is 
in much better agreement with the data. Most modern parton distribution functions, PDFs, 
are evaluated at NLO and fewer new sets are evaluated at LO. However, only LO matrix 
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elements are still employed for some processes such as NRQCD quarkonia production, not 
yet fully calculated at NLO [ |l]]. In addition, the LO matrix elements are inputs to event 
generators like P YTHIA [ ^ . Therefore, LO calculations are still of use, either for speed 
or ease of calculation, if normalized properly relative to NLO. 

Proper normalization is important to avoid over- or underestimating the yields with a 
scaled-up LO calculation. In addition, the shapes of the resulting observable distributions 
may be different because the LO and NLO PDFs are not identical even though they are 
fit to the same data. The LO PDFs must be larger than the NLO PDFs in some regions 
since the NLO total cross sections are generally bigger and the LO PDFs must be larger to 
compensate. In some cases, certain contributions to NLO evaluations are absent at LO. For 
example, y* and vector boson production, denoted by V , proceeds only by qq annihilation, 
q'q V, at LO while at NLO, qq — > Vg and the new process, qg Vq, with a gluon in 
the initial state are also possible. Thus, full NLO calculations, using NLO PDFs, two-loop 
calculations of a, and NLO fragmentation functions, where applicable, are preferable. 

The 'so-called' K factor is used to normalize the LO calculations. The theoretical K 
factor is conventionally defined as the ratio between the NLO and LO cross sections. This 
factor is generally assumed to be constant as a function of the relevant observables. In 
heavy ion physics, the K factor has often been assumed to be 2 without justification. To be 
certain that a LO calculation multiplied hy aK factor does not either over- or underestimate 
the NLO cross section, K should be determined according to a clearly defined prescription 
appropriate to the calculational method. 

There is more than one way to define this factor Before NLO calculations were gen- 
erally available, the LO calculations were scaled up to the data by an arbitrary factor, the 
original K factor. This is an experimental definition, 

^exp = ^. (1) 

If the NLO calculations do not agree with the data, this ratio of data to the higher order 
calculation can also be defined, usually referred to now as 'data/theory', see Ref. [ ^ for a 
recent usage. The K factor is more often determined theoretically. The usual theoretical K 
factor is 

where the superscript (1) refers to the fact that only the first order correction is included. 
When further, higher order, corrections exist, as they do, in part, for Drell-Yan and gauge 
boson production [ as well as heavy quark production [ higher order theoretical K 
factors can also be defined, such as 

„(2) _ CJnNLO „(2') _ CJnnLO 

^th - — ; — ' ^th - ■ 

We do not go beyond K^^^ here. 

There is some ambiguity in the definition of K^^^ . To demonstrate this, we discuss how 
the LO and NLO cross sections are calculated. We first define 

OLo^a(a^) (4) 
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where the exponent n is defined as the power of a, appropriate for the tree-level partonic 
process. We do not include the powers of the electroweak coupling constant here for brevity 
and only concern ourselves with hadroproduction. Thus for Drell-Yan and gauge boson 
production, « = and Olo = o( 1 ) . For direct photon production via qq jg and qg jg, 
n = I and Olo = o(oc.s)- Finally, for heavy quarks produced through the qq QQ and 
gg QQ channels, n—2 and Glo = '^i'^^)- The total NLO cross section is the sum of the 
LO cross section and the next-order correction, o(a"+' ), with one additional power of a^, 

Onlo - Olo + o(a:;+') = o(a';) + o(af ') (5) 

where we denote the total cross section at NLO as Onlo- 

The next order correction for Drell-Yan and gauge boson production, along with virtual 
corrections to qq V, are the new processes qq — > Vg and qg Vq. They are first order 
in Us, so that O^lo = ^(1) +o(ai). The NLO correction to heavy quark hadroproduction 

is a(a^) so that o^lo = ^(f^?) The NLO contribution, a(a"+'), is calculated 

with a, evaluated to two loops and with NLO PDFs. However, there is some ambiguity 
in how to calculate Olo- In principle, at each order of a,, there should be an appropriate 
PDF, LO for a(aj), NLO for a(aj+i), NNLO for a(a^'+2) • • • with a, evaluated to one- 
, two-, and three-loop accuracy respectively. Although possible, at least to NLO, this is 
usually not done. The NLO PDFs and two-loop a, are typically used as defaults for all 
orders of the cross section. In our calculations, we will compare the magnitude and shape 
of the cross sections with Olo evaluated employing LO PDFs and one-loop a^, a full LO 
cross section, and calculations of Glo employing NLO PDFs and two-loop a^. We also 
investigate the effect of these two definitions of Olo on the theoretical K factor for the 
processes considered. 

Thus if Olo is fully LO, we define 

Olo(i) = (6) 
ONLOd) = Oi(a;) + o(a';+i) (7) 

where the subscript '1' refers to the LO PDFs and one-loop a.,. If, instead the LO cross 
section is evaluated with the same PDFs and a.s as the NLO correction, we have 

<7lO(2) = 02«), (8) 

aNLO(2) - 02(al') + a«+>) (9) 

where now the subscript '2' refers to the NLO PDFs and two-loop a,. Recall that the NLO 
correction, a(a"+^), is always calculated with NLO PDFs and the two-loop evaluation of 
a,. 

The total NLO cross section is usually calculated with NLO PDFs and two-loop at 
each order, as in Eq. (^, to determine the size of the next-order correction independent of 
the shape of the PDF and difference in magnitude of ttv. Thus, the most typical way to 
evaluate the NLO theoretical K factor in Eq. (^) is with Eqs. and (H), 

^(1) _ 02(a;)+g(a;+') ^ gNLO(2) 

a2(aj) aLO(2) 
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We have labeled this as the 0* definition because it is most often used by theorists, see 
Ref. [0] for an ( 
Eqs. and (0), 



Ref. [J4|] for an example. It is perhaps more correct to define K^^^ based on cTlo(i)' in 



^(1) _ ai( a;) + o(a;+') _ gNLO(i ) 
oi(a^') ^ aLO(i) 



"'l.Zr ' =^^. (11) 



Now the subscript '1' on Kj^\ refers to the order of the PDFs and a, at which Glo is eval- 
uated. However, if C?nlo(2)' is known but the LO calculation, Olo(i)' is most convenient, 
as in an event generator, it might be more advantageous to define the K factor as 



,(1) _ 02( a;)+o(a;+') _ Onlo(2 ) 
Oi(a?) ^ aLO(i) 



This last definition maximally mixes the PDF and evaluations. Thus ^4,' 2 most de- 
pendent on the difference in the shapes of the LO and NLO PDF fits. Note that, in this case 
only, the K factor cannot be written asK ^ 1 + (tti ). The expressions above are written in 
terms of the total cross sections but are also equally valid for the differential distributions. 

We calculate the differential K^^^, kI^\, and K^^^ with W+, Drell-Yan and QQ pro- 
duction as specific examples. In the case of vector bosons, K^^^ should most strongly 
reflect the difference between the LO and NLO PDF evaluations while QQ production is 
quite sensitive to the order of the ttj evaluation. See Ref. [ |[| for a study of K^^-q in jet 
production. In our calculations, we use the MRST HO (central gluon) [m PDFs in the MS 
scheme for a(a"+^) and Olo(2) ^nd the MRST LO (centi-al gluon) [ ^PDFs for aLO(i) 
unless otherwise specified. 

2. Production in pp Interactions 

We begin our discussion of the K factors with the W+ rapidity distributions. This process 
is independent of the order of since at NLO is always evaluated at two loops. It is 
also the least differential of all the processes considered. We choose the of the three 
gauge bosons because its rapidity distribution is a strong function of the shape of the quark 
PDFs at LO and NLO. We show our results at the top pp energy of the LHC, 14 TeV. The 
dependence of K^^"^ on rapidity is similar for all gauge bosons so that one case is sufficient 
for illustration. See Ref. [ ^ for a comparison of q for all three gauge bosons. We do 
not discuss the pr dependence of either or Drell-Yan production because the pr of 
both are zero at LO. Thus Kj^\pT) cannot be defined for W+ production, only K^^^ [pt]- 

The next-to-leading order, NLO, cross section for production of a vector boson, V, 
with mass my at scale Q^ in a pp interaction is 

^ = H/j Jdxidx2dx5(^^-xxiX2Y(^y-^ln(^^^'j (13) 
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where my is the boson mass, S is the center of mass energy squared, HJ^ is proportional 
to the LO partonic cross section, qiqj V , and the sum over qi runs from u to c. At this 
energy ^ so that the c quark contribution cannot be neglected. The matrices C" and 
C'^ contain information on the coupling of the various quark flavors to boson V. The parton 
densities in the proton are given by fqj{x,Q^) and are evaluated at momentum fraction x 
and scale Q^. For a vector boson without a fixed mass, as in the Drell-Yan process, the mass 
distribution can also be calculated by adding a dM in the denominator of the left-hand side 
of Eq. (jTsj) and the delta function 5(M — my) on the right-hand side. 

The prefactors //^^ are rather simple [ ^]. In the case of W+ production. 



2k Gf 
3 V2 S 



(14) 



where Gf = 1.16639 x lO^^ GeV^ and mw — 80.41 GeV. For virtual photon production 
via the Drell-Yan process, there are three contributions to Hjf: virtual photon exchange, Z" 
exchange, and y* — Z*^ interference. 



Ht 



f-zP 



H. 



Hi 



9M25' 

a Gf m\ (1-4 sin^ 9^ ) (M^ - m|) 
9^~S (M2-m|)2+m|r| 
1 Gf m2 otzF; 



e,|(l-4|e,-|sin2ew) 



3^/2 5 (M2-m|)2+m|r| 



(l + (l-4|e,|sin2 9w)^) 



(15) 



where sin^Gw = 1 -m^/m|, niz = 91.187 GeV, Tz = 2.495 GeV, and Tz^i+i- = 2>361%. 
Note that here we use / = j since / and j must be the same flavor for neutral gauge boson 
production. Now Hf^ depends on the pair mass and goes inside the mass integral. 

The functions A,j(x, Q^) in Eq. ( pj| ) are universal for all V [ We work in the MS 
scheme. The NLO correction to the q'q channel includes the contributions from soft and 
virtual gluons as well as hard gluons from the process q'q Vg. We have, up to NLO [ H], 
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iNLO 



(^,e') 



-x) 



(16) 
(17) 
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»NLO 
99 



, 1- 



-4(l+x)ln( % ) -8(l+x)ln(l-x) 



(18) 



1 -X 

1 



lnx + 5(l -x) 



1 



In 



6 In 



16 



ln(l 



8C(2)-16 



l-x 



The LO contribution is a delta function. The NLO contribution is proportional to ttj, cal- 
culated to two loops with rif = 5 active flavors. The first three terms of A^i'*^(x,2^) are the 
real contributions from qq Vg while the last three terms are the soft and virtual gluon 
contributions from the qTI vertex correction. The general integral of the 'plus' functions 
from soft gluon emission in Eq. (|l6|) is [ |lO| 



dxf{x) 



ln'(l -x) 



1 



= dx 



/(I 



ln'(l-x) + :!^ln'+^(l-fl) 



1 



The quark-gluon contribution only appears at NLO through the real correction qg 
At this order [ E[, 



A«(x,e^) 



a.(e2) 



(19) 

^qV. 
(20) 



871 



2(l+2x^-2x)ln 



(1 



\2n2 



l-lx^ + 6x 



Note that, since at NLO, this contribution is only real, no delta functions or plus distri- 
butions appear. Using the delta functions in Eq. (|l3|) we find xi 2 — \J 2^/x5'exp(±3') at 
NLO. For A^°(x, 2^) and the terms in A^°(x, 2^) contributing to the vertex correction in 



all terms proportional to \n{Q /my) 



Eq. (|l|), x'j 2 = ^/Q^/Sexp{±y). If we take ^ 1 
drop out. 

We now define the coupling matrices, C" and C'^, in Eq. (pj|). The superscripts rep- 
resent the initial (i) and final (f) state quarks or antiquarks while the arguments indicate 
the orientation of the quark line to which the boson is coupled [ Q]. The couplings 
are elements of the CKM matrix. They are nonzero for 0^{qk,qi) if + e; = ±1 and for 
0^{qk,qi) if = ±1 + e/. In both cases, they take the values jV^j^, Following Ref. [ ^], 
we take Vi,d — cos 9c ~ Vcj and Vus ~ sin 9c ~ — Krf with sin 9c ~ 0.22. 

The q'q convolution in a pp interaction is 



I _f,P,ixuQ%'^ix2,Q^)C"{qi,qj) 

- cos2 9c te(xi , e2)/|(x2, e^) , e')/.^(x2, e^) 



(21) 



sm^dc[fl!ixi,Q^)ffix2,Q^)+f^{xi,Q^)fP{x2,Q^)]+[xi ^X2 



For the qg channel, we have 

E _(/.^(^l,G')/;(-^2,e') + [xi^X2] ]C\q„qt) 



(22) 
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+ [xi ^ X2] . 

The coupHngs do not enter explicitly in the qg convolution because each distribution is 
multiplied by (cos^ 9c + sin^ 9c). 

In Fig. |l] we show the W+ rapidity distributions calculated with Eqs. (||)-(^. The 




Fig. 1. The rapidity distribution in pp collisions at 14 TeV evaluated at Q = mw+ ■ The 
solid curve is cyNLo(2)^ the dashed curve is Olo(i)' the dot-dashed curve is Onlo(1) ^nd the 
dotted curve is 0lo(2) ■ 

'standard' NLO calculation, C7nlo(2)' is shown in the solid curve. The accompanying LO 
result, CJlo(2)' is given in the dotted curve. The two curves are nearly parallel to each other 
over most of the available rapidity space. Both increase slowly with rapidity until y ~ 3 
where the turnover point is reached. The NLO calculation increases slightly faster than 
the LO, most likely due to the qg contribution, /,f (xi)/|'(x2), since the gluon distribution 
increases with decreasing x when x is small. The increase with y is due to the /,f (xi )/j (^2) 
component of the qq annihilation cross section. The valence u distribution in the proton, 
flivi^i)^ is increasing until xi ^ 0.126, corresponding to the peak of fuyixi) at y ~ 3. 
Likewise /-^{xi) is rising with decreasing X2- Only when fuyixi) is falling with growing 
xi, beyond the peak, does the cross section begin to decrease. 

The same effect is seen, albeit more strongly, for C?nlo(i) ^i^d ^LO(i)- We see that 
Olo(2) increases 6% for < y < 3 while C5lo(i) increases by 18% over the same interval. 
The difference arises primarily because f-j{x) increases faster with decreasing x at LO than 
at NLO. Note also that C?lo(2) is about 20% larger than (Jlo(i) at y = 0. This difference 
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is due to the fact that the NLO PDFs are larger in the region around x ~ 5.6 x 10^^. As 
y increases, the LO valence and sea distributions become closer to the NLO distributions, 
so that, in the tail, aLo(i) > C'lo(2)- This trend is also reflected in the corresponding NLO 
distributions. 

The resulting K factors are shown in Fig. ^. The 'standard' K factor, ^t/Jq, is the 




Fig. 2. The three K factors for W+ production. The solid curve is K^^ q, the dashed is K^^ \, 
and the dot-dashed is k[^\. 

smallest, ~ 1.13, for < y < 3.5 where it begins to increase. The rise is due to the greater 
importance of the fg{xi)fg{x2) contribution as X2 0. The factor Kj^^\ is only slightly 
higher, ^ 1.16, and also rather flat over most of the rapidity range. The largest K factor 
is K^^\, as might be expected from inspection of the curves in Fig. [l|, due to the greater 
difference between Onlo(2) ^"d Olo(i)- The decrease for y > 2 is due to the steeper rise of 
the LO rapidity distribution as the peak of fuy (x) is approached. The minimum of A^'J 2 is 
reached when Onlo(2) drops below Onlo(i)- 

We have also checked the dependence of K on scale and PDF. The scale dependence 
enters not only in the PDFs and a, but also in the logarithms 111(0^ /niy) in Eqs. ( [l6| ) 
and (po|). When = niy, as in our calculations so far, the three scale dependent terms 
in A^{x,Q^) drop out and A^^'^{x,Q^) only depends on ln{{l - x)^ /x). If > m^, 
the logs are positive, enhancing a(as), but if < niy, the scale-dependent logs change 
sign, decreasing o(aj). Scale evolution increases the low x density of the sea quarks and 
gluons while depleting the high x component and reducing the valence distributions. The 
X values do not change when is varied. Thus the higher scales also tend to enhance the 
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cross sections. At these large values of Q^, does not vary strongly so that increasing or 
decreasing by a factor of 16 is not a large effect. The total cross sections are larger when 
— Amy and smaller when = my/ A, opposite the scale dependence at lower energies, 
due to the larger parton densities at low x and high Q^. The same trend was also seen in 
Ref. [ ^]. We find that all three defintions of the K factor exhibit the same dependence. 
Increasing gives a 7% larger K factor while decreasing reduces the K factors by 
~ 10% for the MRST set. 

The variation of the results with PDF is also significant. The MRST and CTEQ [ 
pT] ] results are similar but the CTEQ cross sections are a few percent higher. In this case, 
the Olo(2) ^nd Onlo(2) distributions are not as broad as Olo(i) ^"d Onlo(i)- The peak of 
aLo(2) is at y < 3 while aLO(i) peaks aiy ^ 3.2. Because of this difference, ^4! 9 is not 
independent of rapidity anywhere but decreases almost to unity where Olo(i) ^ ^nlo(2)- 
When the GRV 98 distributions [ [ij] are used, the rise with rapidity is much stronger than 
that in Fig. |l} However, now the differences between Olo(i) ^lo(2) rather small, 
aLo(i) peaks about 0.25 units of rapidity higher than Olo(2) and with a slightly larger value. 

In this case, the K factors are all smaller than for MRST, k[1\ - k[1\ - 1.09. Now even 
^lii 2 °"ly ^ percent different than the other two definitions over two units of rapidity. 

3. Drell-Yan Production in Pb+Pb Interactions 

We now discuss dilepton production in the Drell-Yan process where a, also does not en- 
ter until NLO, as in the previous section, but the LO cross sections are more differential, 
depending on pair mass, M, and rapidity. Since, at collider energies, the Drell-Yan contribu- 
tion to the dilepton continuum is small compared to heavy quark decays and will probably 
not be cleanly separated, we discuss the Drell-Yan yield at the CERN SPS. The Drell-Yan 
cross section is of interest here because it is the only source of lepton pairs in the mass 
region above the \(/' for \/5 =17.3 GeV. As such, it has been used to normalize the 
yield and determine its suppression as a function of transverse energy, Ej- (See Ref. [ |l3| 
for more details.) The NA50 collaboration [ |lj] measures the dilepton continuum above 4 
GeV in PbH-Pb collisions and fits the experimental K factor to the dilepton yield in the mass 
range 4 < M < 9 GeV by a LO Drell-Yan calculation with the MRS A [ ^ NLO PDFs. The 
Drell-Yan contribution below the //\|/ peak, in the range 2.9 < M < 4.5 GeV for NA50, is 
then calculated at LO and multiplied by this same experimental K factor, assuming that it 
is independent of mass and rapidity in the interval Q <y < \. It is appropriate to check this 
assumption. 

The Drell-Yan cross section as a function of mass and rapidity is also given by Eq. ( pj| ) 
with the delta function added for the mass, as explained earlier. We do not include any 
nuclear effects on the PDFs. The convolution of the prefactors , Eq. (p3|), with the 
parton distribution functions in a nuclear collision in the qq channel is: 

£ ('/^(xi,e2)^A^(_,2,e2)^^A.(^^^g2)^A.(_^^ g2A (23) 
qieQ V / 
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= < ({Z^/,f (xi,e2) +NAf:{x,,Q^)} {Zsffix2,Q^) +NBmx2,Q^)} 
+ 2ABfP{xuQ^)f^{x2.Q^)) 

+ Hf, (^{ZAfSixuQ') +NAf^{xuQ')} {zBf^{x2,Q') +NBf^{x2,Q^)} 
+ 2ABfP{xi,Q^)fP{x2,Q^)) + [xi ^X2,A^B] . 



We let i/j^ represent all charge +2/3 quarks and H^^ all charge —1/3 quarks. In the qg 
channel the convolution is 



We have taken the isospin of the nuclei into account, denoting the proton and neutron 
numbers in a nucleus with mass A by Za and Na respectively. 

Since we are interested in the makeup of the dilepton continuum in heavy ion collisions 
at the SPS, the values of \/S and M relevant for our discussion are both considerably lower 
than in the previous section. We first present the Drell-Yan mass distributions in Fig. ^ 
At \/S =17.3 GeV, the y* — Z and Z" contributions to the Drell-Yan cross section, though 
included, are negligible. Note also that the slope of Olo(i) with M is different than Olo(2)- 
In particular, aLO(i) is similar to and larger than 0lo(2) for M > 4 GeV, where NA50 fits 
their data to the LO calculation, while, for M < my, aLo(2) is larger than cyLo(i)- Thus, 
choosing a NLO PDF to fit an experimental K factor to a LO calculation could lead to some 
differences in the calculated Drell-Yan cross section below the J/\\t mass region. 

The calculated K factors, shown in Fig. ^ as a function of M, demonstrate different de- 
pendencies on M. At this lower energy, all K factors are larger than for production in 
the previous section. Some of this difference is because = <C so that a, is larger 

for low mass Drell-Yan production. The standard K factor, Aj'Jg, increases from ^ 1.75 

to 2.3 between 1 and 10 GeV while K^^\ is nearly constant at ^ 1.95 and Kj^^2 decreases 
over the entire mass interval. None of these values agree with the constant experimental K 
factor of 1.7 [ |l4|| integrated over M obtained by NA50 with the MRS A [ |l5|] set. Since 
the Drell-Yan data have relatively low statistics, a precise determination of the Drell-Yan 
cross section is difficult for M > 4 GeV, especially when the Pb-nPb data are divided into 
centrality bins. Thus, any measure of changes in the slope of the mass distribution in Fig. || 
is unlikely in nuclear collisions at this energy. However, as shown in Figs. |] and ^ the 
relative slopes and thus the theoretical K factors do change. A constant K factor cannot 




(24) 



+ {ZBf^ix2,Q^) +NBmx2,Q^)}+2AfnxuQ^)) 
+ Hf, {{ZAf^ix, , +A^^/«(xi , Q^)} 

+ [zBfl{x2,Q^) +A^B/^(x2,e2)} +2A/f(x,,e2)) 



+ [xi <-> X2,A «-> B 
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2.5 5.0 7.5 10.0 



M (GeV) 

Fig. 3. The Drell-Yan mass distribution in Pb+Pb collisions at 17.3 GeV evaluated at 
Q — m. The solid curve is cTnlo(2)' the dashed curve is cyLo(i)' 'he dot-dashed curve is 
Onlo(I) ™d the dotted curve is cTlo(2)- 

account for this. The ratios of the to Drell-Yan cross sections used to represent the 
suppression pattern as a function of Ej, whether in the measured Drell-Yan or in the mini- 
mum bias analyses, both depend on extrapolation of Olo(2) multiplied by the experimental 
K factor 

The results for Olo(2) ^nd Onlo(2) have also been calculated with the MRSA low 
PDF [ ITsI] used by NA50 [pi]. Since there is no corresponding LO PDF for the MRSA set, 

(1) 

we only compare K^^^ ,j. We find that Onlo(2) with MRSA is smaller than that with MRST 
for all masses, as much as 23% smaller at M = 10 GeV. On the other hand, aLO(2) is a 
few percent larger with MRSA for M < 6 GeV, dropping below the MRST cross section 
at higher masses by up to 16%. These differences are relatively small but, taken together, 
reduce k[1\ for MRSA to ^ 1 .8 at 4 GeV, closer to the experimental K factor of NA50. 

The Drell-Yan rapidity distributions and K factors are shown in Figs. ^ and ^ respec- 
tively. We have separated the results into three mass bins: 2.9 < M < 4.5, the region in 
which NA50 compares the Drell-Yan and data; 4.5 < M < 6 GeV; and 6 < M < 9 
GeV, above the region. Of the latter two bins, the lower mass bin should correspond 
to the best statistics for NA50. Unfortunately, the NA50 Drell-Yan data are not binned in 
rapidity and are only given in the interval Q <y < 1 . Note that the rapidity range becomes 
more restrictive as the mass increases. Only the lower two bins reach y =\. 

The trends of the cross sections as a function of mass in Fig. || are shown more clearly 
in Fig. |. In the lower mass bin, Onlo(2) > Onlo(i) and aLO(2) > Clo(i) for y < 0.75. In 
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Fig. 4. The three K factors for Drell-Yan production at the SPS as a function of mass. The 
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the intermediate mass bin, the two ways of calculating the LO and NLO cross sections 
coincide up to y w 0.5 where the results calculated with (Tlo(i) become larger. Finally, in 
the large mass bin, Onlo(i) > <5nlo(2) and c?lo(i) > c?lo(2) for the entire rapidity range. 
The K factors as a function of rapidity reflect the cross section pattern. Again the K factors 

are aU greater than 1 .7 and, except for l^^, are not constant but tend to grow with rapidity, 
especially K'f^^y 

Interestingly, when the cross sections are calculated with MRSA, the corresponding 
A'j^'q are always somewhat smaller than those with MRST and ~ 1.7 near j w in 
the lowest mass bin. The agreement with NA50's experimental K factor is more likely a 
fortuitous accident than a precise fit. Given the rather limited Drell-Yan statistics of NA50 
for M > 4 GeV, it is doubtful that a better measure of the experimental K factor could have 
been obtained. It would be worthwhile, however, to calculate the full NLO cross section, 
either (Tnlo(i) or (Jnlo(2) depending on the chosen PDF, and obtain the experimental K 
factor at NLO. Our results suggest that it would be significantly smaller but stiU above 
unity, ~ 1.2. 

4. Heavy Quark Production in Interactions 

We now turn to evaluations of the heavy quark cross sections and their K factors. We 
concentrate only on the bare distributions rather than introduce empirical fragmentation 
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functions which, for heavy quarks, remain at a rather primitive level [ |l^. We can expect 
the biggest variations in the K factors here because the LO cross section is akeady pro- 
portional to a^. The order at which is calculated is important because the lower Q^, 
proportional to = nig + pj in distributions, means is much larger for charm and 
bottom quarks than for gauge bosons. Differences in Aqcd between LO and NLO PDF 
fits result in similar values of at one and two loops when the LO and NLO values of 
Aqcd respectively are used. However, large discrepancies may result if the LO value of 
Aqcd is used in a two-loop evaluation of a.,. Thus, one must be careful to use PDFs and ttj 
evaluations that are compatible. 

Heavy quarks are produced at LO by gg fusion and qq annihilation while at NLO g{q + 
q) scattering is also included. To any order, the partonic cross section may be expressed 

ik I) 

in terms of dimensionless scaling functions f-j ' that depend only on the variable r\ = 
V4m2 -1 [|, 

a,jis,mlQ') = ^ f: {4:zaM))'i (n) In' ( ^] , (25) 

'"e k=o 1=0 \ '"s / 

where s is the partonic center of mass energy squared, niQ is the heavy quark mass and 
is the scale. The cross section is calculated as an expansion in powers of a, with 
k = corresponding to Olo = C^lct.?)- The first correction, k ~ I, corresponds to o(oCs). 
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Note that no distinction is made between LO and higher evaluations of ttj in Eq. (|25|). 
The scale is generally assumed to be the same for both the renormalization, and 
factorization, Q^, scales since all PDF analyses make the assumption that — Qp = Q^. 
It is only when k> 1 that the dependence on Q\ can be distinguished since when k — 1 
and I = 1, the logarithm \n{Q^ /m^) appears. The dependence on appears already at 
LO in the argument of in the partonic cross section and in the parton densities. The 
NNLO corrections to next-to-next-to-leading logarithm with k = 2 have been calculated 
near threshold [ |[|. The complete calculation only exists to NLO. The total hadronic cross 
section is obtained by convoluting the total partonic cross section with the PDFs, 

appiS,ml) = jU t —f.'{xi,Q^)ff (-,Q^) Oij{zS,ml,Q^) (26) 

where the sums / and j are over all massless partons. 

To illustrate the dependence of the total QQ cross sections upon mn and Q^, we show 
the cc and bb cross sections as a function of mass and scale in Figs. 0and ||. The charm 
mass is varied between 1 .2 GeV and 1.8 GeV in a/S = 200 GeV pp collisions in Fig. |7[ We 
choose this energy because it is close to that of the PHENIX charm measurement at RHIC 
[ |l7|], -y/S =130 GeV. A more complete measurement should be coming from the 200 GeV 
data taken in the latest run. The bottom mass is varied between 4.25 GeV and 5 GeV for 
pp colUsions at \/5 = 41.6 GeV, the energy of the HERA-B experiment which has recently 
presented a measurement of the bb total cross section [ [l^]. Unlike charm production at 
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Fig. 7. The cc total cross sections as a function of charm quark mass mpp collisions at 200 
GeV. The lefthand plot shows the results for = while the righthand plot is calculated 
with = 4m^. The solid curves are Onlo(2)' th^ dashed curves Olo(i)' '^he dot-dashed 
curves cyNLO(i) ^"d the dotted curves cyLO(2)- 



RHIC, bottom production at this energy is in the near-threshold region where resummation 
techniques can be applied [ ^ . 

The cc cross section as a function of mass in Fig. |7|is calculated for two different values 
of the scale, = rrt^ (left-hand side) and 4m^ (right-hand side). Note that only scales 
greater than are shown because = m^. /4 < 1 Ge for all charm masses considered, 
below the minimum scale of most PDFs. Thus results for scales lower than would not 
be particularly meaningful. For this particular energy, the scale dependence is not strong, 
10-20% at LO for = 1 .5 GeV, and 40% at NLO. The stronger scale dependence at NLO 
is unusual and shows that charm could be difficult to treat at higher orders because of its 
relatively low mass [ |l9|]. The mass dependence is stronger for the NLO calculations than 
the LO and is also stronger for the larger scale. Note that in both cases displayed in Fig. ^ 

C^NLO(2) < C^NLO(i) because cTlo(i) > (5lo(2)- The difference is small since A'j'J' > 2 for 
all three definitions of the K factor Thus CT(aJ) alone is significantly larger than either 
calculation of the LO cross section. The most important sources of the difference in the 
LO cross sections are the larger LO gluon density and the dominance of the gg production 
channel. The charm K factors tend to be larger for lower masses and smaller scales. 

We have shown the results for one specific energy here. However, the scale dependence 
changes with energy, as already noted in Section^. The growth of the cc total cross sections 
with energy is slower for = m^. At lower energies, o{m^) > a{4m^). As the energy 
increases, a(4m^) > o{m^). The value of \/S at which this cross over occurs changes with 
Mc, VSfulOO GeV for Mc = 1 .2 GeV and 3 TeV for ot,. = 1 .8 GeV. The change in 

energy dependence of the cross section with scale is due to the evolution of the PDFs at 
large and low x, similar to what was observed for production in Section ^. 

In Fig. ^ the bb cross section is given as a function of mass for three different values of 
Q^, ml/4- (upper left), (upper right) and 4m| (bottom center). The larger bottom quark 
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Fig. 8. The bb total cross sections as a function of bottom quark mass in pp collisions at 
41 .6 GeV. Clockwise from the upper left, the results are calculated with = mj^, and 
4ml respectively. The soUd curves are (Tnlo(2)' the dashed curves (Tlo(i)> the dot-dashed 
curves crNLo(i) the dotted curves C7lo(2) • 



mass now allows us to calculate the cross section for scales lower than m| and still obtain 
meaningful results. In this case, the scale dependence is reduced at NLO relative to LO. 
The mass dependence is stronger at NLO, as for charm, and increases somewhat with scale. 
The heavier quark mass reduces the impact of C7(a^) on the total cross sections. Indeed, 

for = C?lo(1) > ^nlo(2)- This effect is probably due to the one-loop evaluation 

in <Tlo(i) since the one-loop a, is larger than the two-loop a.,, particularly at low scales. 

The bottom K factors are all smaller than those for charm, K^^^ < 2 in most cases. Thus bb 
production, at least at this energy, is better under control than charm. Like charm, however, 
the K factors also are larger for lower masses and higher scales. The scale dependence 
of bb production with energy is smaller than for charm. The bottom cross sections follow 
the hierarchy a(m|/4) > C7(m|) > a{4ml) for all energies and orders. Near threshold the 
difference in cross sections due to scale can be nearly a factor of 10 but by \/S =14 TeV, 
the difference is only a few percent. The smaller scale dependence is also an indication that 
bottom production is better under control than charm. 
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The mass and scale parameters used in our further calculations are determined by 
obtaining a 'best' fit to the data without an experimental K factor Using the MRST PDFs 
to calculate Onlo(2)^ we found =1.2 GeV and — Arn^ and mt, — 4.75 GeV with 
= for bottom [ ^ |2|]. Although there is a great deal of scattered data up to 

= 63 GeV on the cc total cross section, our choice for the bottom parameters are in line 
with conventional wisdom rather than data. When other PDFs were used, the parameters 
favored varied somewhat [ |TJ |2|] but, in the following, we compare results from different 
PDFs using the same values of mg and Q^. 

Since the charm mass obtained from these evaluations is somewhat smaller than the 1 .5 
GeV mass generally used and the NLO calculations with mi, = 4.75 GeV underestimate the 
Tevatron pp bb data [ ^sjl, it is reasonable to expect that higher order corrections beyond 
NLO could be large. Indeed, the HERA-B cross section [ |l8|] agrees with the NNLO-NNLL 
cross section in Ref. [ |[|, suggesting that the next order correction could be nearly a factor 
of two. Thus the NNLO correction could be nearly as large as the total NLO cross section. 

Unfortunately, the NNLO-NNLL calculation is valid only near threshold. The pp data 
at higher energies, while not total cross sections, also show a large discrepancy between 
the perturbative NLO result and the data, nearly a factor of three [ ^3|]. This difference 
could be accounted for using unintegrated parton densities [ although these unintegrated 
distributions vary widely. Another, more mundane explanation to the bb discrepancy is an 
incomplete understanding of the hadronization process [ ^ . If some resummation is needed 
in the high energy regime to account for the Tevatron data, there remains the difficulty of 
connecting the regimes where near-threshold corrections are applicable and where high- 
energy, small X ^ niQ/^/S physics dominates. This problem is increased for charm where, 
even at low energies, we are far away from threshold. 

Our method is perhaps the most straightforward-using Onlo(2) ^nd ignoring higher- 
order corrections to fit the data. This is not difficult for cc because the data are extensive. 
However, there are less bb data. The "KT p cc data tend to favor lighter charm masses. It 
is difficult to say if the same is true for bb. A value of nib = 4.75 GeV, which underpredicts 
the Tevatron results compared to NLO cross sections [ ^3|], agrees reasonably well with 
the average of the iiTp data. However, for the HERA-B measurement to be compatible 
with a NLO evaluation, the b quark mass would have to be reduced to 4.25 GeV, a value 
which might be more compatible with the Tevatron results. Therefore, a full NNLO calcu- 
lation with ntc ^1.5 GeV and mi, ~ 4.75 GeV might agree with the QQ data without an 
experimental K factor. A more quantitative statement is not possible. 

We now employ our inferred values of mg and to calculate the K factors as func- 
tions of energy and PDF. Figure ^ illustrates the danger inherent in calculating, Olo(i) ™d 

then multiplying either by o °^ ™ arbitrary factor of 2: the cross sections can be con- 
siderably overestimated. Indeed, K factors are generally applied to avoid underestimating 
the cross sections. We see that A'j'J q is generally greater than 2 at low energies, drops some- 
what and finally increases with energy. This dependence of the standard K factor could be 
attributed to more inherent theoretical uncertainties in low and high energies where dif- 
ferent resummation techniques are applicable. There is a somewhat stronger effect for bb 
smce near the threshold region, until a/S 100 GeV. On the other hand. 
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Fig. 9. The three K factors for QQ production in pp colUsions as a function of energy. The 
cc results, calculated with nic = 1.2 GeV and Qp' = are given on the lefthand side, 
while the bb, calculated with = m\ and mi, = 4.75 GeV, are on the righthand side. From 
top to bottom the results are given for the MRST, CTEQ5, and GRV 98 PDFs. The solid 
curves are q, the dashed are and the dot-dashed are 2- 



charm production is well above threshold for all energies shown since 2mc/VS < 0.1 for 
\/S > 20 GeV. Note that kI^\ and K^^^^ both decrease with energy with a larger decrease 
for cc (lefthand plots) than bb (righthand plots). This decrease can clearly be attributed 
to decreasing x ^ Q^/S with increasing V5. The X range of charm production is from 
jc ~ 0.1 at 15 GeV to lO"* at 14 TeV. The LO gluon PDF is larger than the NLO gluon PDF, 
particularly at low scales, so that Olo(i) increases faster than C7nlo(2) ^th 2 decreases 
to near or below unity at high energies. The values of x and ^ are larger for bb production 
so that the decrease in /T^' 2 with energy is slower, k'^^ does not decrease as quickly as 

k[^\, and is almost independent at high energies, but is still smaller than K^f^. Thus at 
high energies, 0lo(i) <^nlo(2) in some cases and including a K factor on Olo(i) would 
overestimate the cross section by a factor of 2 — 3. The same trends are observed for results 
with the CTEQ5 and GRV 98 PDFs. 
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The PDF dependence is illustrated in Fig. [T^ with the ratios of the CTEQ5 and GRV 98 
LO and NLO cross sections to the MRST cross sections. The ratios depend most strongly 




10^ lo3 10"* 102 lo3 10'*^ 
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Fig. 10. Ratios of the total QQ cross sections calculated with the CTEQ5 (top) and GRV 
98 (bottom) PDFs to those calculated with the MRST PDFs as a function of energy. The cc 
results, calculated with = 1.2 GeV and = 4m^, are given on the lefthand side, while 
the bb, calculated with = and mt, — 4.75 GeV, are on the righthand side. The solid 
ratios are Onlo(2)> the dashed ratios C7lo(i)' the dot-dashed ratios C?nlo(i) ^"d the dotted 
ratios Olo(2)- 

on the relative gluon PDFs. The CTEQ5 and MRST NLO gluon distributions are very 
similar so that the ratios of Olo(2) ™d CTnlo(2) ^re nearly unity without much variation 
over the entire energy range. There are slight differences for bb at low energies due to the 
larger qq contribution near threshold. The GRV 98 NLO gluon distribution is less similar 
to the MRST gluon, leading to a larger variation of the ratios with energy. The ratios 
involving the LO gluon PDFs are larger and more energy dependent because the LO gluon 
distribution is not as well constrained. The biggest differences are in the ratios of aLO(i) 
since the large gg contribution to o(oc^) reduces the effect of the LO gluon distribution in 
the Onlo(i) ratios. 

To complete our discussion of heavy quarks, we now turn to the charm and bottom dis- 
tributions and their corresponding K factors, shown in Figs. |ri|jl4| . The necessity of NLO 
evaluations of QQ production is clear if one wants to study pair production and correlations. 
The pT of the QQ pair is zero at LO and only becomes finite at NLO. LO calculations are 
thus only useful for obtaining the single quark distributions in pj, rapidity and Feynman 
X, Xf, and for the QQ pair, the y, xp and invariant mass, M, distributions. Therefore the 
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Fig. 11. The charm quark pj (a), rapidity (b), and cc pair mass (c) for pp colHsions at 
200 GeV. The soHd curve is <7nlo(2)' the dashed curve is CTlo(1)' the dot-dashed curve is 
*^NLO(i) ™d the dotted curve is cyLo(2)- 



quantities for which we calculate the K factor, besides the total cross sections, are the sin- 
gle quark pr and rapidity distributions and the pair mass distribution. We will evaluate 
the differential K factors for charm distributions at RHIC, ^/S = 200 GeV, and bottom at 
the LHC, V5 — 5.5 TeV, both appropriate energies for ion colliders. We calculate Olo(i) 
analytically, resulting in smooth distributions. However, C7lo(2) ^nd o(oc^) are calculated 
with the NLO Monte Carlo code described in Ref. [ 25 1 which uses a two-loop evaluation 
of tti as a default. The finite statistics of the Monte Carlo result in some fluctuations in the 
distributions, particularly apparent in the K factors. 

The charm calculations at \/5 = 200 GeV for RHIC correspond to the total cross sec- 
tions shown in F ig. y\ with ==1.2 GeV and — 4m^. It is clear that the LO pT dis- 
tributions in Fig. |1 l|have a steeper slope than the NLO distributions. The charm rapidity 
and pair mass distributions, on the other hand, are more similar. The K factors increase 
rapidly with pj in Fig. even though we have used m|. as a more appropriate scale for 
the distributions [ This increase can be expected since ©(ocj) is inherently at higher 
Pt due to contributions from gg — > gg* gQQ where the QQ pair is opposite a gluon jet. 
Contributions such as these begin to be important when pj > mg, as is obvious in Fig. |2[ 
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Fig. 12. The three K factors as a function of charm quark pj (a), rapidity (b), and cc pair 
mass (c) for pp collisions at 200 GeV. The solid curve is g, the dashed is k'^\, and the 
dot-dashed is /ci'l. 



The increase in the K factor with pr is also only important when x is relatively small — the 
K factors for bb production at RHIC increase very slowly with pr over the same interval. 
The pT dependence of the K factors clearly should be taken into account when scaling up 
a LO calculation. Conversely, the K factors are nearly constant with rapidity and pair mass 
which are pj- -integrated and are weighted by the low pj part of the cross section. The same 
trend with rapidity seen in Figs. ^ and || is also observed here: ^4,' i) increases slowly with 

rapidity while k[^\ decreases. The exact behavior at large rapidity is difficult to determine 
because the Monte Carlo fluctuations are largest near the edges of phase space. 

A careful examination of the K factors may reveal some slight difference in magnitude 
between Figs, [l^and^ at the same energy. This shift is due to the somewhat different scales 
used to calculate total cross sections and distributions. For total cross sections, the only 
relevant scale is <^ nig in the ln(2^/mg) term in Eq. (p5|). However, when calculating 
distributions, pj- -dependent logarithms also enter which need to be controlled by the mj 
scale. Shifting the scale from °^ ntg to °^ nij has the effect of modifying the total cross 
section obtained by integrating the distributions. This shift in scale is manifested by a 
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Fig. 13. The bottom quark pj (a), rapidity (b), and bb pair mass (c) for pp collisions at 5.5 
TeV. The solid curve is cyNLO(2)' the dashed curve is cTlo(i)' the dot-dashed curve is Onlo(i) 
and the dotted curve is cTlo(2)- 



change in the magnitude of the K factors as well. 

The results for bb production in a/S = 5.5 TeV pp collisions at the LHC are shown in 
Figs. 13 and 14 The K factors clearly begin to increase for pr > 5 GeV where pr > mi,. 
The increase with pj is slower than for charm production. The K factors are essentially 
constant over the rapidity range shown since the edge of phase space has not yet been 
reached aty — 5. Also, away from threshold, 2mi,, the K factors are constant as a function 
of mass. Note that K^^^q exhibits the strongest threshold effect near M ^ 10 GeV. 



5. Summary 

We have discussed three ways of defining the theoretical K factor to NLO. The first, Kj^\, 
is most appropriate for determining the theoretical uncertainties between cross sections of 
different orders because it allows the most straightforward determination of the next-order 
effects. In cases where the use of the LO cross section is necessary for speed of calculation, 
such as in event generators, it is most appropriate to make a full LO evaluation and then 
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Fig. 14. The three K factors as a function of bottom quark pj (a), rapidity (b), and bb pair 
mass (c) for pp colHsions at 5.5 TeV. The soHd curve is k[^^q, the dashed is and the 



dot-dashed is K, 
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th,2- 



multiply by either k'^^^ or K^^2- 

We have shown that the theoretical K factor is only approximately constant away from 
threshold and away from the edges of phase space. Thus for high energy W+ production, 
K is only independent of rapidity for y <2. The Drell-Yan K factor is not constant with 
mass, increasing at low masses and also changing at high masses, increasing or decreas- 
ing depending on which definition is used. The rapidity range over which K is constant 
for Drell-Yan production at the SPS decreases with increasing mass as the edge of phase 
space is approached at lower rapidities. The K factors for heavy quarks are most strongly 
dependent on pj since the NLO corrections are large for pj > niq. Similarly, a strong 

dependence o^K^^^q on pj in inclusive jet production was found in Ref. 

Clearly a constant K factor is inappropriate for all kinematic variables. To best utilize 
knowledge of the next-to-leading and higher order corrections, the differential K factor 
should be determined as a function of the kinematic variables of interest, such as pj, y or 
M. Whether k[^-^ or ^4^' 2 ^^^^ with C7lo(i) is somewhat a matter of taste. However, it 
should perhaps be kept in mind that the PDFs are evaluated to NLO using Onlo(2) ■ Thus 
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2 perhaps the most relevant for extending full LO calculations to NLO. Therefore 
there is no clearcut answer to the question posed in the title. Whatever definition is chosen, 
it is obviously important that the determination of Ka^ be clearly described and applied 
consistently in a given calculation. 
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ful discussions. 
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